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AN EFFICIENT METHOD FOR THE SOLUTION OF 
SCHWINGER DYSON EQUATIONS FOR PROPAGATORS. 

MARC P. BELLON 



Abstract. Efficient computation methods are devised for the perturbative 
solution of Schwinger-Dyson equations for propagators. I show how a simple 
CL(, computation allows to obtain the dominant contribution in the sum of many 

^ ' parts of previous computations. This allows for an easy study of the asymp- 

totic behavior of the perturbative series. In the cases of the four-dimensional 



^\ . supersymmetric Wess-Zumino model and the (^5 complex scalar field, the 

^Sl ' singularities of the Borel transform for both positive and negative values of 

the parameter are obtained and compared. 

1. Introduction 

D \ In preceding works [H 12]) we have shown how to solve Schwinger-Dyson 

f^, ' equations for propagators. This allows to compute the /3-function in models 

without vertex divergencies, e.g., the supersymmetric Wess-Zumino model and 

CN ' some version of a six-dimensional theory of a scalar field with cubic interactions. 

The resolution method in the first of these works [1] is exact, but does not 

Q"^ \ allow for an easy understanding of the properties of the computed perturbative 

series. This method is also computationally heavy and would scale badly for 
Schwinger-Dyson equations with a greater number of propagators. 

^O I The approximate differential equations introduced in [2] have the double 

advantage of allowing an easy derivations of the asymptotic behavior of the 
perturbative series and of a complexity which is only quadratic in the pertur- 
bative order. However, as can be seen in the second equation proposed for 
- the Wess-Zumino model, this simplicity is lost if we have to get more precise 

rN ' results. 

C^ ■ In this work, I propose a procedure which allows to systematically improve 

the second method, while remaining computationally simple and making the 
asymptotic behavior of the perturbative series transparent. In particular, I shall 
derive the dominant singularity of the Borel transform of the perturbative series 
on the positive axis, through the added precision on the asymptotic behavior 
of the perturbative series. In this work, I limit myself to the contributions of 
the dominant poles, which are sufficient for the asymptotic behavior, but the 
inclusion of additional terms would allow to reach any desired precision, for a 
fraction of the computational cost of the method used in p.J. 

All these computations are based on two techniques, which have been de- 
scribed extensively in []4[2j, elaborating on the work of Kreimer and Yeats [5]: 
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a renormalization group equation allows to deduce the full propagator from the 
renormalization group functions and a Mellin transform of the one loop diagram 
gives the renormalization group functions from the propagators. These results 
will be freely used here. 

In the present paper, we will see in a first step how the contribution of a simple 
pole in the Mellin transform of the diagram can be computed recursively, allow- 
ing the computation of an infinite sum of derivatives of the propagator while 
sidestepping the computation of the individual derivatives. This is immediately 
applied in the following section to the solution of the linear Schwinger-Dyson 
equation for a scalar model in dimension 6 which was first solved in [3]. The 
new method surely recovers the same numerical values and is computationally 
comparable, but it is far superior for the asymptotic analysis of the perturbative 
series and it prepares for the more complex examples. 

In the cases of non- linear Schwinger-Dyson equations, some terms cannot 
be computed with this first method and I extend it to the case where the de- 
nominator depends on the sum of some variables of the Mellin transform. This 
method is then applied to the supersymmetric Wess-Zumino model. The main 
result thus obtained is a precise asymptotic study of the perturbative series: 
beyond the dominant contribution which gives a pole for the Borel transform 
of the series, the exact nature of the singularity on the positive axis, which was 
discovered in the numerical computation of [1], is uncovered. Finally, the scalar 
model in dimension 6 is studied in the case of a non-linear Schwinger-Dyson 
equation. Again, a singularity on the positive real axis is predicted, but it would 
have been difficult to spot numerically, since the non alternating part is smaller 
by a factor bigger than n'^. Finally, I conclude on the possible extensions of 
this work, in particular with respect to the inclusion of higher loop corrections 
in the Schwinger-Dyson equations. The singularities of the Mellin transforms 
for higher loop diagrams are poles with denominators of the type envisioned 
here: the dominant contributions of the high order derivatives of the Mellin 
transform can be computed by the methods presented here. 



2. Contribution from a simple pole. 

We keep the fundamental insight in [2j that the Taylor series is dominated 
by the contribution of the singularities near the origin. In the one-loop case, 
the Mellin transform is expressible through F functions and these singularities 
are single poles. The first case to consider is therefore the one of a simple pole 
l/{k + x). 

The corresponding contribution to the 7- function is: 






p=0 ^ ^ p=0 ^ ^ 

In this expression, 7^ is the p derivative with respect to the logarithm of the 
impulsion L of the correction to the propagator G{L). However, the seemingly 
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local sum of derivatives in equation ([T|) can be given an equivalent integral form: 
(2) Ik= I G{L)e'''^dL 

J —oo 

If k is negative, the integral has to be taken from to +00. For /c = 1, we 
recover the integral 

/ G{p')d{p% 
Jo 

which appeared in the work of Broadhurst and Kreimer [3]. 

The point is that I can be efficiently evaluated through the use of the renor- 

malization group. The action of 7 + /3da on G{L) gives its derivative with 

respect to L; 



(3) 



1-0 

(7 + ^ada)Ik = / dLG{L)e''^dL = G(0) - kik 
J —00 



The last equality comes from an integration by part. Equation ([3|) can be used 
to obtain an efficient recursive computation of the perturbative expansion of I. 

3. The linear case. 

As a first application, let us go back to the linear Schwinger-Dyson equa- 
tion for a 6-dimensional scalar field first studied in [1] . The one- variable Mellin 
transform is a rational function and a simple solution method was presented 
in [4J. This amounts to acting on the Schwinger-Dyson equation with an ade- 
quate differential operator in the logarithm of the momentum L to reduce the 
right-hand side to a constant. In |2j, I proposed to act with a simpler dif- 
ferential operator. This allowed in particular to more easily account for the 
asymptotic behavior of the perturbative solution, but the truncated right-hand 
side precluded an exact solution. 

With the insight of the preceding section, it is possible to obtain the exact 
solution in a form which allows for an easy asymptotic analysis. Indeed, through 
the use of I2 and Z3 which represent the effect of the propagator on the single 
poles 1/(2 + x) and 1/(3 + x), equation (10) of [2J can be written: 

(4) 7 + 7(2aa„-l)7 = a(Z2-X3). 

The recursive evaluation of X2 and X3 through the use of equation (J3]) , completed 
by the identity /3 = 27, allows for a rapid evaluation of the perturbative series 
for 7. The result coincides with the one stemming from the partial differential 
equation of [1], providing a check of the computation. 

It is then easy to obtain from equation (JH the asymptotic behavior of the 
perturbative series. Let us fix the notation for the perturbative coefficients of 
7, Z2 and I3: 

00 00 CXD 

(5) 7 = ^c„a^ X2 = J^d„a", l3 = Y,fna''. 

n=l n=0 ra=0 

It is convenient to convert the term bilinear in 7 to the form {ada — 1)7^, so 
that with the first values ci = 1/6 and C2 = —11/6^, equation ^ gives the 
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following relation: 

2n 22(n - 1) , , 

(6) C„+i ~ a'^''^ 63 Cn-l+dn-fn- 

In the equation ^ , the dominant components come from the cases where either 
7 or Ifc is of maximal order, so that one obtains: 

(7) dn ^ --c„ ——dn-l 

/o^ f ~ _i 2n- 1 ^ 

\°) In ~ qC^ Jn— 1 

Now, it is consistent to suppose that dn and /„ are proportional to c„ and replace 
Cn-i in the right hand side of the previous equations by the its approximate 
value —-Cn to solve for the unknown factor. One obtains: 

2 D 

Putting all together, one obtains the following asymptotic relation: 

(10) c„^., - - ^3 + ^J Cn, 

which is well verified on the exact solution. By considering additional terms in 
the expansion of 7^ and the recursive definition of dn and /„, one could obtain 
additional terms in the expansion in - of the ratio of c„ over c^-i- 

4. General poles. 

When studying non-linear Schwinger-Dyson equations, Mellin transforms 
with many variables must be used. However, the obtained singularities are 
really simple, since they correspond to the divergences appearing when some 
subgraph becomes scale invariant: the singularities are simple poles whose de- 
nominators only depend on the sum of the Mellin variables associated to the 
subgraph. 

The formulation of the pairing between propagator and Mellin transform 
appears quite dissymmetric, but it is in fact totally symmetric. The evaluation 
of a function of the propagator G{L) can be reduced to the following pairing of 
G{L) with the Mellin transform H{x): 

00 ^ 

(11) Y.-^LG{L)\L=od:H{x)U=o 

^ — ' n! 

n=0 

For simple enough Mellin transforms H, it will be more convenient therefore to 
write this duality in reverse order, i.e., as some differential operator acting on 
the propagator G{L). The structure will simply be: 

(12) ^(9l)G(L), 

since for an analytic H, the Taylor series gives the value of the function. The 
situation becomes more interesting when there are many variables, since a sum 
of derivatives with respect to Li can be converted to a simple derivative if we 
identify the different variables. The methods of section [2] can then be applied 
to evaluate the effect of a simple pole. 
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However, one must pay special attention to the numerator: if in one variable, 
any numerator can be reduced to a scalar, it is not the case in a multidimensional 
setting. It is even possible, if we start from a multiloop diagram, that the 
numerator itself be a meromorphic function with poles. The identification of 
the different variables Lj has therefore to be done after the application of the 
differential operator corresponding to the part of the numerator which depends 
on the variables in the denominator. The final result is that we can evaluate 
easily the following type of pairings: 

1 " 

(13) Z= N{dL,....,dL:)Y^G{L,\ 

^ Li=-=L„=0 

We obtain an equation similar to eq. ([3]), apart that the anomalous dimension 7 
is multiplied by the number of fields n and that G(0) is replaced by an expression 
involving A^. 



(14) [k + n-i + Pada)l = N{dL„. • • , 5l J J] G{L, 

5. WeSS-ZuMINO MODEL. 



=i„=0 



Equipped with these new tools, let us consider the super symmetric Wess- 
Zumino model. It is now possible to compute additional terms of the asymptotic 
behavior of the perturbative series. In particular, I shall obtain the exact nature 
of the singularity on the positive real axis. The aim is to include all the poles 
at unit distance from the origin, that is to use the approximation of the Mellin 
transform introduced in equation (21) of [1|, corrected to have the exact xy 
term: 

(15) H.. ,) = (1 + .,) (-1- + -1- - 1) + 1—2^ + i.,. 

To solve the Schwinger-Dyson equation with this approximate Mellin trans- 
form, we need two functions in addition to 7 itself: F, the sum of the 7„ asso- 
ciated to 1/(1 -|- x) and H, the function associated to the term xy/{l — x — y). 
For the same approximation, the differential equation obtained in section 3.4 
of [2] is very complex. 

The renormalization group yields to the following functional equations for F 
and H: 

(16) F = 1 - 7(3aaa + 1)F, 

(17) H = 72 + 7(3aa„ + 2)F. 

Using the approximate Mellin transform (jlSp . the Schwinger-Dyson equation 
then takes the form: 

(18) 7 = 2aF -a- 2a'y{F - 1) + ^a{H - 7^) 

The three equations (|16ll7|18p can be converted in recursion equations for the 
perturbative expansions of the functions 7, F and G. 
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The numerical solution is easy. Compared with the complete calculations 
of [l], the obtained precision is good, with errors which are always less then 1 
percent and are asymptotically of the order of half a percent. This divides the 
error by a factor 10 with respect to the cruder approximation of [2j and more- 
over capture the non-alternating part in the development of 7: an asymptotic 
analysis of this system of equations will allow to obtain the exact nature of the 
singularity on the positive axis in the Borel transform of the series for 7. 

Let us write 'j = Yl ^noP'-, F = Yl fnO-^ and H = Y hnO^- It is easy to show 
that ci = 1, C2 = —2, /2 = 1. With the rapid growth of all the coefficients, one 
obtains the following asymptotic forms of the equations ()16|17p : 

(19) fn+i ^ -(3n + l)/„ + 2(3n-2)/„_i-c„+i, 

(20) /i„+i ~ 2c„-4c„_i + (3n + 2)/i„-2(3n-l)/i„_i + 8cn_i. 

In the first of these equations, (3n — 2)/„_i and c„+i are respectively propor- 
tional to — /„ and 2/„, yielding the asymptotic relation: 

(21) /n+i~-(3n + 5)/n. 

Remarking that the dominant term in Cn is 2/„_i, one recovers the dominant 
asymptotic relation for c„, Cn+i — — (3n + 2)c„, observed in [1] and proved 
in [2]. For /i„, one obtains similarly the asymptotic relation 

(22) hn+i =^ 3n/i„. 

This term is non-alternating and is the dominant contribution to the singularity 
of the Borel transform on the positive real axis. The ratio of the growth of fn 
and gn is in absolute value 1 — 5/3n. This translates in a ratio n~^'^ between the 
absolute values of these two contributions: this ratio coincides nicely with the 
numeric results and confirm the subleading character of the c„ term in eq. (120p . 



6. The 0g model. 

In this case, one considers the action of di + 9£ on the Schwinger-Dyson 
equation. The left-hand side is then 

(23) 7 + 7{3ada - 1)7, 

and the right-hand side is based on the following Mellin transform: 

r(2 + x)r(2 + y)r(l-x-y) 



(24) Hix,y) 



r(l-x)r(l-y)r(4 + x + y) 



I approximate H with the contributions from the poles for x or y equal to —2 or 
—3 and the one for 1 — x — y = 0. For y = —2, the residue is (1 — x)(l — x/2) and 
is converted to {l+xy/2){l-\-xy/4:) not to mar the low order terms. For y = —3, 
the residue has two more factors and is — (1 + x){l — x){l — x/2)(l — x/3). For 
l — x — y = 0, the residue is xy/12(l + xy/2). We therefore will use the following 
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approximation of H: 

^ '^^ ^ ^ 2 ^^ 4 ^ \2 + x 2 + y 2 

/. xy.,^ xy.,^ xy.. xy. f 1 1 1 
-(1 ^)(1 + — )(1 + — )(1 + — ) \ 

(25) +-^X2/(1 + ^)- i ^xy. 

12 2 1 — X — y 21b 

The last term in this equation adjust the coefficient of xy to the exact value, 
so that the cubic term in 7 is exact. Most remarkably, this coefficient appears 
as the sum of the contributions from poles more distant from the origin, with 
poles in A; + x giving negative contributions and the poles in k — x — y giving 
positive contributions. For the terms in x'^y or xy'^ where both contributions 
add up, the remainder is larger. 

For the asymptotic analysis, these additional terms do not matter. We need 
three functions in addition to the anomalous dimension 7, Z2 and I3 as in the 
case of the linear Schwinger-Dyson equation and a function H associated to 
the term xy{l + xy/2)/{l — x — y). We have similar equations for X2 and Z3 
than in the linear case: 

(26) 2Z2 = l-7(l + 3aaa)X2, 

(27) 3Z3 = l-7(l + 3aaa)X3. 

The recursion equation for H just has additional source terms with respect to 
the Wess-Zumino case, due to the more complex numerator: 

(28) H = ^^ + \-tl + 7(2 + ^ada)H. 

Finally, using the approximate form of the Mellin transform (|25p and neglect- 
ing terms divisible by (xy)^ which are subdominant, one obtains the following 
equation: 

7 + -ii^ada - 1)7 = a(2X2 -\) + a|7(2X2 - 1 - ^7) + «h2(2^2 - I - H 

-a(2X3 - i) - ^07(2X3 -l-h) 

(29) +^«72(2X3 - I - §7) + h^H - ^a^\ 

This gives the first two non-zero coefficients of 7, ci = g, C2 = —j^- The 
dominant contribution for c„+i comes from the second term of the left-hand 
side and is proportional to n/2. The dominant terms for X2 and X3 are in fact 
exactly the same than in the case of the linear Schwinger-Dyson equation, since 
the changes from 2n to 3n in the recursions for c„ and in equations (J26l27p have 
compensating effects. Eq. ([9]) remains valid at the dominant level. At the next 
to leading order, we have therefore: 

n 13 

(30) Cn+i ~ -(3n + l)(cic„ + C2C„__i) -h 2dn - 2/„ ~ -(- -F — )c„. 

The recursion operation for H psp will differ from the Wess-Zumino case 
through the values of the first coefficients of 7. The source terms are not 
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important for the non-alternating component of H. We obtain 

(31) hn+i ^ -(3n + 2)hn - j^i^ri - l)/i„_i ~ {^ - Y^)^"- 

The absolute values of the factors in the recursions for /i„ and Cn are in the 
ratio 1 — 31/(9n), meaning that the absolute value of /i„ is smaller by a power 
31/9 = 3,444. . . of n than Cn- Compounded with the fact that /i„ intervenes 
only in c„+i, this non-alternating component in c„ could easily have been missed 
in a numerical study. In a sense, the situation is even worse, since the 7^ term 
in equation (j28p will give a larger contribution, so that hn itself is dominated 
by this alternating component. However, this is not really a problem since 
the asymptotic recurrence relations are linear and the full solution will be a 
superposition of the particular solution proportional to c„_i and a solution of 
the equation ([3T]) . 



7. Conclusion. 

In this work, the solution of Schwinger-Dyson equations for propagators has 
made a new step forward: I reach not only the leading asymptotic behavior 
of the perturbative terms, but the subleading, "wrong sign" contribution, as 
well as systematic corrections in powers of 1/n. These computations can be 
made easily more precise through the inclusion of additional poles of the Mellin 
transform or the adjunction of a few monomials: a suitable combination of 
the two methods should allow to obtain the solution of the Schwinger-Dyson 
equations with very high precision for a fraction of the computational cost of 
the methods of [T]. However, such a full precision computation is not really 
useful since the Schwinger-Dyson equations I consider here are but the first 
approximation to the full system of Schwinger-Dyson equations. 

In the computation of higher order corrections to the Schwinger-Dyson equa- 
tions, the power of the present methods should be precious. Indeed, in a higher 
loop primitively divergent diagram, the number of individual propagators is 
higher and the full evaluation through the straightforward methods of [Ij, much 
more complex. For a Mellin transform in k variables, the number of different 
derivatives up to total order n scales as n^~^, with the evaluation of each of 
these terms at perturbative order 2n again implying approximatively n^~^ op- 
erations. Some time for space bartering could reduce somewhat this growth by 
reusing partial products, but in any case, this complexity would make any com- 
putation, whether explicit or asymptotic, very complex, even if the evaluation 
of the derivatives were not difficult per se. 

In the approach presented here, the situation is much more manageable: the 
pole structure of the Mellin transform, which is linked to the divergences of sub- 
diagrams, is of the form analyzed in section H] and their contributions can be 
computed in linear time. Furthermore, one can identify the residues which are 
important for the asymptotic behavior of the perturbative series and focus the 
analytical evaluations on them. This will be the subject of a forthcoming pub- 
lication [3|. If, up to now, these works have dealt only with Schwinger-Dyson 
equations for propagators, the efficiency of the methods introduced here should 
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be important in the more challenging cases of systems including Schwinger- 
Dyson equations for vertices. 

In the spirit of p], we could also deduce systems of differential equations 
for the anomalous dimension and some auxiliary functions, with the hope of 
obtaining information on the asymptotic behavior of the anomalous dimension 
at large coupling. However, the singularity of the Borel transform on the pos- 
itive axis indicates that the Borel resummation is not uniquely defined. It is 
therefore a challenge to determine whether such a system of differential equa- 
tions, which would be determined from purely perturbative considerations, can 
be given a non-perturbative meaning. 
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